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Electron-electron (e-e) and electron-hole (e-h) interactions are often associated with many exotic 
phenomena in correlated electron systems. Here, we report an observation of anomalous excitons at 
3.75 , 4.67 and 6.11 eV at 4.2 I< in bulk- SrTiOa. Fully supported by ab initio GW Bethe-Salpeter 
equation calculations, these excitons are due to surprisingly strong e-h and e-e interactions with 
different characters: 4.67 and 6.11 eV are resonant excitons and 3.75 eV is a bound Wannier-lifce 
exciton with an unexpectedly higher level of delocalization. Measurements and calculations on 
SrTii_ x Nba ,03 for 0.0001<a;<0.005 further show that energy and spectral-weight of the excitonic 
peaks vary as a function of electron doping ( x ) and temperature, which are attributed to screening 
effects. Our results show the importance of e-h and e-e interactions yielding to anomalous excitons 
and thus bring out a new fundamental perspective in SrTiC> 3 . 


I. INTRODUCTION 

Strong electronic correlations in transition metal ox¬ 
ides are known to drive rich phenomena such as su¬ 
perconductivity, colossal magnetoresistance and metal- 
insulator transitions [1]. An example of a model tran¬ 
sition metal oxide is the perovskite-type SrTiOa (STO). 
Indeed, STO in various bulk forms and recently its het¬ 
erostructures with other oxides, have shown exotic phe¬ 
nomena such as superconductivity [2-6], magnetisms [5— 
9], metal-insulator transitions and two-dimensional elec¬ 
tron gas [10-14]. However the role of electronic cor¬ 
relations, particularly electron-hole (e-h) and electron- 
electron (e-e) interactions, remains unclear. In the past, 
it was assumed that the role of e-h and e-e interactions 
in STO was negligible, if any, and thus most calcula¬ 
tions were done without implicitly incorporating such 
many-body effects [15-18]. Recently, a theoretical study 
by solving the ab initio Bethe-Salpeter equation (BSE) 
within the GW approximation has shown that if such 
many-body effects are present in STO, this might yield 
to the so-called resonant excitonic effects and a predic¬ 
tion of a resonant exciton at ~6.4 eV, which then should 
all be manifested in its complex dielectric function within 
a broad energy range [19]. If this is true, it makes STO 
even more interesting and important to study. Until we 
solve this fundamental problem, we may not be able to 
understand, for instance, the origin of such exotic phe¬ 
nomena as well as to resolve discrepancies in previous 
optical spectroscopic studies particularly near absorption 
edge characteristics in STO [20-27]. An ideal tool to 


study resonant excitonic effects is spectroscopic cllipsom- 
etry (SE) because it probes neutral excitations, e.g. e-h 
pairs and excitons, providing important information in 
particular on two-particle correlations. Indeed, SE has 
been used and was successful to unravel the presence of 
resonant excitonic effects in graphene [28, 29] as theoret¬ 
ically predicted [30, 31]. 

In this work, systematic temperature-dependent SE 
measurements supported by theoretical ab initio BSE 
calculations on top of a GW electronic structure (GW- 
BSE) are used to unravel the presence of anomalous ex¬ 
citonic effects and thus to identify the role of e-h and e-e 
interactions in SrTii_ x Nb x 03 for 0<x<0.005. Measure¬ 
ments and calculations on SrTi 1 _ a ,Nb x 03 are performed 
to determine effects of charge carrier screening on the 
optical spectra. 

II. SAMPLE PREPARATIONS, 
EXPERIMENTAL METHODS AND 
THEORETICAL CALCULATIONS 

SrTi0 3 (100) and SrTi 1 _ x Nb x 0 3 (100), x = 0.0001, 
0.0005, 0.001 and 0.005 samples procured from Crystec 
are used for the measurements. All samples are single 
side polished of size 10 mmx 10 mm x 0.5 mm. Atomic 
force microscopy (AFM) measurements show that the 
rms roughness for all the samples are less than 5 A. 

Spectroscopic ellipsometry (SE) is a technique that 
measures the polarization state change of light upon re¬ 
flection from (or transmission through) a sample. The 
measured ellipsometric angles ’L and A are related to 
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the material properties of the sample by the complex re¬ 
flectance ratio 


q = r p /r s = tan’Pe lA (1) 

where r p and r s are the Fresnel reflection coefficients 
for p and s polarized light, respectively [32]. 

Spectroscopic ellipsometry is performed with a com¬ 
mercial rotating analyser ellipsometer (J.A. Woollam, 
Inc. V-VASE) covering the energy range of 0.6 - 6.5 eV. 
Low temperature measurements are performed within a 
cryostat (Janis) and at base pressures in the 10~ 9 torr 
regime. Liquid helium and liquid nitrogen are used in 
the open cycle cryostat. The angle of incidence for all 
cryostat measurements is 70°. Major focus is put on 
the pseudo-dielectric function (PDF) as it can be di¬ 
rectly obtained from the measured ellipsometric quan¬ 
tities (dq A) and it is identical to the complex dielectric 
function (£ 1 , 62 ) for an isotropic bulk sample with a per- 
fecly smooth top surface [32]. With an optical model 
assuming a perfectly flat substrate having infinite thick¬ 
ness, the PDF ((ei), (£ 2 )) can be obtained from the mea¬ 
sured (’P, A) directly. In the experimental set up SrTiC >3 
(STO) and SrTii_ x Nb x 03 samples are measured in air 
and vacuum and hence the optical model is taken as an 
air/sample or vacuum/sample interface. In this case 
the PDF can be calculated using the following expres¬ 
sion [32]: 


(s) = sin 2 0 i 


1 + tan 2 0 i 



( 2 ) 


Where 9i is the incident angle and g is related to the 
measured (\P, A) by the expression (1). 

The GW-BSE calculations are performed for x = 0 
and 0.005 emphasizing the role of e-h and e-e interac¬ 
tions, particularly on near edge spectra, as well as screen¬ 
ing effects in electron doping case. Note that neither a 
Kramers-Kronig transformation nor optical modelling is 
required here to further ensure that the observed opti¬ 
cal phenomena are truly intrinsic property of the system 
under study. 

To solve the static BSE problem, it can be convenient 
to adopt a basis of independent-particle (IP) transitions 
[33] . In this basis the excitonic eigenvalue problem reads 


= E x A\ (3) 


FIG. 1. (Color online) Experimental critical points and 
comparison with theory for STO. (a) (ei) and (£2} for the 
lowest (4.2 K) and highest (350 K) measured temperatures, 
(b), (c) Temperature dependent differences 5x{£i) = (ei)(T) — 
( £i)(T = 350A'),i = 1,2. (d) Plots of £1 and £2 calculated 
using GW-BSE and GW-RPA (all graphs are blue-shifted by 
47 meV to match experimental data). The green vertical 
dashed line are determined using second derivative analysis. 


where t and t! label the IP-transitions, A\ is the eigenvec¬ 
tor relative to the eigenvalue E\ . The BSE absorption 
spectrum is determined calculating the macroscopic di¬ 
electric function: 



E\ — u> — irj 


( 4 ) 
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FIG. 2. (Color online) GW-BSE analysis of the peak 
Exi. Peak Exi: E\ = 3.644 eV. GW-BSE spectrum and 
IP-spectrum are represented by the red and green lines re¬ 
spectively. The histogram of is presented in blue. The 

cumulant weight function is reported in dotted black. A black 
arrow indicates the energy E\ = 3.644eV. 


FIG. 3. (Color online) GW-BSE analysis of the peak Ex 2 . 
Peak Exa: E\ = 4.607 eV. GW- BSE spectrum and IP- 
spectrum are represented by the red and green lines respec¬ 
tively. The histogram of is presented in blue. The cu¬ 

mulant weight function is reported in dotted black. A black 
arrow indicates the energy E\ = 4.607eV. 


where Vo(q) is the Coulomb interaction, 77 is a positive 
infinitesimal quantity and, 

Pvc(q) = We- Iqr |c> (5) 


where c and v are conduction and valence states respec¬ 
tively. The peaks are found for exciton energies E\ (the 
zeroes of the denominator in Equation (4)). The weight of 
the excitation at Ex is given by the squared modulus of 
a sum over IP-transitions (the numerator of in Equation 
(4)), namely 


rroo 


E^(q) 


with Sl(q) = pt(ai) A \ 


( 6 ) 


where S\ (q) are the IP-weights and the sum is made over 
all transitions of the system (converged spectrum), which 
is stressed by the 00 symbol. 

For a given peak at energy E\ , one can define an IP- 
transition energy interval I = [E IP — S, E IP + 5] and 
look at the partial spectral weight I J2tei ^(q)!- This 
gives information about the absolute contribution to the 
spectral weight of all transitions of energy close to E IP . 
By choosing many energy intervals, a histogram can be 
obtained. 

However, the IP-weights <S^(q) are complex numbers 
which can sum in constructive or destructive way depend¬ 
ing on their phase. In order to account for destructive 
contributions, the partial spectral weight 


(q) = 


t:t<E IP 


Si(q) 


(7) 


is defined as the spectral weight obtained summing only 
IP-transitions of energy lower or equal to E IP . With the 


help of this quantity, the cumulant weight 


fx(E IP ) 


lim 

q —>0 



can be plotted as a function of E IP . 


( 8 ) 


III. RESULT AND DISCUSSION 


The temperature dependence reveals a conspicuous 
evolution of the PDFs. In Fig. 1(a), for the sake of 
clarity, the PDFs of only the highest (350 K) and low¬ 
est (4.2 K) measured temperatures for STO are com¬ 
pared where considerable sharpening of features is ob¬ 
served at the lowest temperature. The most important 
observations are revealed in the temperature difference 
data 6 T (£i) = ( £i)(T ) - ( £i )(T = 350 K),i = 1,2 in Figs. 
1(b) and 1(c). These results reveal features around 3.35, 
3.80, 4.12, 4.70, 4.82, 5.40 and 6.24 eV at 4.2 K (Fig. 
1). The structure at 3.35 eV corresponds to the indirect 
band-gap, which is in agreement with previous reports 
[20, 22, 27], 

The main observation is two sharp peaks at ^3.80 
and ^4.70 eV as clearly seen in <5 t(£2 ) (Fig. 1(c)). In- 
triguingly, they both are the most sensitive to temper¬ 
ature, with decreasing temperature their intensities in¬ 
crease monotonically. Furthermore, St{£ 1 ) shows a sharp 
change in phase at these structures (Fig. 1(b)) due to 
Kramers-Kronig relations between £\ and £ 2 - This fur¬ 
ther supports that these peaks are intrinsic optical prop¬ 
erties of STO. Detail fitting based on the second deriva¬ 
tive of dielectric function reveals that the these particular 
peaks are of purely excitonic nature with energy positions 
at 3.75 and 4.67 eV at 4.2 K (see discussion below). For 
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FIG. 4. (Color online) GW-BSE based analysis of excitons Exi and Ex 2 . Excitonic wavefunction 4/(rh, r e ) for (a) Exi and 
(b) Ex 2 on the ac Ti-O plane (red spheres: oxygen, cyan spheres: titanium). The probability |T(rj,, r e )| to End at position 
r e the electron belonging to the exciton when its bound hole (black sphere) is on the central oxygen atom is indicated by the 
colour scheme from blue to red. 


the rest of discussions, we will refer to these excitons as 
Exi and Ex 2 , respectively. 

We perform two different theoretical calculations, 
random-phase approximation (RPA) on top of GW (GW- 
RPA) and GW-BSE calculations, to determine the STO 
complex dielectric response £\ + ie 2 (Fig. 1(d)). Note 
that while the GW-RPA calculations do not account for 
e-h interactions, the GW-BSE implicitly incorporates e-h 
and e-e interactions. Intriguingly, the GW-BSE calcula¬ 
tions are in remarkable agreement with the experimen¬ 
tal data, especially the presence of structures at 3.75 eV 
and 4.67 eV, which are absent in GW-RPA calculations. 
This directly implies that the Ex! and Ex 2 features cor¬ 
respond to distinct peaks identified as excitons in the 
GW-BSE spectrum. Furthermore, comparison of GW- 
RPA and GW-BSE results reveal a significant spectral 
weight transfer, a fingerprint of correlations, and hence 
the crucial role of excitonic effects to explain experimen¬ 
tal spectra of STO. More thorough discussions on theo¬ 
retical calculations are given below. 

DFT and G 0 W 0 calculations have been performed 
with the free plane wave package ABINIT [34], Norm- 
conserving pseudopotentials including semicore electrons 
have been employed for Sr (10 electrons in valence) 
and Ti (12 electrons in valence). DFT calculations, 
performed within the local density approximation, con¬ 
verge with a basis set cutoff of 70 Ha. The density 
has been computed by sampling the Brillouin zone with 
a llxllxll Monkhorst-Pack fc-point grid, while the 
Kohn-Sham eigenvalues have been computed on a grid 
of arbitrarily shifted 512 fc-points (8x8x8). The G 0 W 0 
self-energy has been computed using a cutoff energy of 
25 Ha (2400 plane waves) for both the exchange and the 


correlation terms. The later has been constructed on 
two frequencies following the Godby-Needs plasmon pole 
method where 200 bands have been included in the sum 
over states. The G 0 W 0 corrections have been evaluated 
on 35 fc-points (8x8x8 Gamma-centered mesh) and their 
values have been interpolated to obtain the correction on 
the shifted grid. 

Optical spectra (IP and GW-BSE) have been obtained 
using Go Wo energies as input. The dielectric matrix has 
been computed on the 512 shifted fc-points grid, with a 
cutoff energy of 40 eV (33 plane waves) and including 40 
bands. In GW-BSE calculations, the screened Coulomb 
interaction has been computed with the same parame¬ 
ters used for the correlation self-energy. Diagonalization 
of the excitonic Hamiltonian has been performed either 
in the iterative Haydock scheme with 150 iterations, or 
in the full diagonalization if also the eigenvalues where 
needed for analysis. 

In the case of SrTii-^Nb^Oa, x = 0.005, in order to 
investigate the evidence of increase of screening due to 
doping, we solve GW-BSE using G 0 W 0 energies of the 
undoped case but accounting for the doping only in the 
screened Coulomb interaction. This has been computed 
with the same parameters as the undoped case, but start¬ 
ing from LDA eigenvalues and eigenfunctions obtained 
including an extra charge of 0.005 electrons per unit cell. 

Figures 2 and 3 show how the two excitonic peaks (Exi 
and Ex 2 ) are qualitatively different. It has to be noted 
that the theoretical e 2 obtaned using GW-BSE calcula¬ 
tions has been plotted with a global blue shift of 47 meV 
for comparison with the experimental result in Fig. 1(d) 
of the main text. Here in the following discussions we 
refer to the same corresponding peaks but with the asso- 
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FIG. 5. (Color online) Double-derivative data and fit of (a) 
(ei) and (b) (£ 2 ) for STO at 300 K. The fitting parameters 
can be found in Table I. 


dated theoretical energy positions. The theoretical peak 
positions are 3.644 and 4.607 eV for Ex! and Ex 2 re¬ 
spectively. In the case of the peak Exi at 3.644 eV (the 
onset), only IP-transitions of higher energy are summed. 
This can be clearly seen in Fig. 2 by looking at the his¬ 
togram of Sf (blue impulses): it is different from zero 
only starting from energies E\ > 3.87 eV. This allows 
for identifying this exciton as a bound exciton since its 
binding energy is well defined. Moreover, the first con¬ 
tributions are also the most important in absolute value. 
The analysis shows that the cumulant spectral weight 
(black dotted line) immediately increases as soon as the 
histogram has non vanishing contributions. On the other 
hand it decreases gradually as long as the absolute contri¬ 
butions become smaller: therefore all IP-transitions sum 
in a constructive way, and all participate to the spectral 
weight of the first exciton. However, at higher energy 
(between 6 and 8 eV) the cumulant function goes above 
1, having a maximum for 7 eV. This means that if we in¬ 
clude IP-transitions up to 7 eV, we will overestimate the 
intensity of the first peak. Its weight is then reported 
to the converged result, by including higher energy IP- 
transitions, which sum with an opposite phase, reducing 
the spectral weight and taking the cumulant function to 
the correct asymptotic value of 1. By comparing the cu- 
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FIG. 6. (Color online) Experimental critical points in the 
PDFs of STO and SrTii-^NbcOs. (a) (ei) and (e 2 ) for STO 
and SrTii-^Nb^Os with different doping levels, (b), (c) Dop¬ 
ing dependent differences Sd{£i) = (£i)(x) — {si){x =0),* = 
1 , 2 . 


mulant function with the IP-spectrum (green line in Fig. 
2), one can see which structures of the IP-spectrum are 
mixed to create the exciton at 3.644 eV. In this case, 
contributions are mostly from structures between 4 and 
6 eV resulting in a relatively strong excitation effect. 


The analysis of the peak Ex 2 leads to a different inter¬ 
pretation. As it can be seen in Fig. 3 by the histogram 
S* (in blue), IP-transitions of energy E\ < 4.644 eV are 
included in the summation, which makes it difficult to 
give a clear-cut definition of its binding energy. How¬ 
ever, the fact that the cumulant weight function (dotted 
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TABLE I. Parameters of the double derivative fit of (e) for STO at 300 K. 


Critical Points 

E (eV) 

A 

T (meV) 

<t> (deg) 

Type 

E 0 

3.350(29) 

0.27(1) 

156(31) 

303.62(38) 

2D 

Exi 

3.801(4) 

0.43(3) 

193(4) 

86.45(6) 

Excitonic 

E x 

4.147(16) 

4.63(85) 

277(19) 

311.96(15) 

2D 

EX2 

4.712(19) 

8.40(2.34) 

531(21) 

174.24(5) 

Excitonic 

E2 

4.860(1) 

8.39(2.95) 

264(20) 

195.06(3) 

2D 

e 3 

6.000(338) 

16.63(6.76) 

1537(356) 

73.49(31) 

2D 

Ex 3 

6.113(10) 

0.58(8) 

325(14) 

255.62(10) 

Excitonic 


black) is vanishing, tells us that these transitions have op¬ 
posite phase and they sum destructively. It is interesting 
to note also that the most important absolute contribu¬ 
tions (that is the highest blue impulses) have opposite 
phases. This is evident from the small peak in the cumu- 
lant function: it raises because of the first impulse, but it 
is immediately killed by the transitions of slightly higher 
energy (the second impulse) which have opposite phase. 
Amongst the transitions summing constructively in the 
peak Ex 2 , the bump at 7 eV (green line) contributes for 
~ 20 %. 

Based on GW-BSE calculations, the origin of the var¬ 
ious structures in the optical spectra of STO can now 
be interpreted more comprehensively. Analysis of the 
exciton composition in terms of the remixing of single¬ 
particle e-h transitions confirms the strong excitonic na¬ 
ture of both peaks (see above). The dominant contribu¬ 
tions to Exi are coming from O -2p—> Ti -3d t 2g transi¬ 
tions in the range from the band-gap to ^4.5 eV. Thus, 
Exj can clearly be identified as a bound exciton. The 
calculated binding energy for Exj is ~220 meV, which 
is in good agreement with our experimental estimation. 
In contrast, even though Ex 2 shares the same O -2p^ 
Ti -3d 1 2 g orbital character as Exi, the dominant contri¬ 
butions are from larger energies (from ~ 4.5 to ~ 5 eV). 
Therefore, Ex 2 can be interpreted as a resonant exciton. 

The real space plot of the excitonic wavefunction T(r^, 
r e )| as calculated by BSE for the bound exciton Exj is 
shown in Fig. 4(a). The shape of the exciton wave- 
function clearly has the characteristic of the Ti -3d t 2g 
electrons, consistent with the above-mentioned compo¬ 
sitional analysis. However, surprisingly the bound exci¬ 
ton Exi appears as a highly delocalized Wannier-lzfce ex¬ 
citon emerging from strongly localized Ti -3d electrons. 
Meanwhile, the resonant Ex 2 exhibits a more localized 
character than the bound Exi, which is contrary to the 
conventional Wannier picture [35], as seen in Fig. 4(b). 

Finally, a feature labelled Ex 3 is observed at 6.11 
eV (this energy position is determined using double¬ 
derivative analysis as shown below. Based on BSE cal¬ 
culations at higher energy bands [19], this feature has 
been attributed to resonant exciton via 0-2p—> Ti -3d e g 
transitions with spectral weight mostly from states at 


higher energy of 8-10 eV. In order to analyze the detailed 
temperature dependence, it requires optical spectra with 
much higher energy, which is beyond our current mea¬ 
surement. Nevertheless, this further supports the impor¬ 
tance of electronic correlations in STO. Furthermore, the 
prominent features observed experimentally at 4.12 eV, 
4.82 eV, 5.40 eV are inter-band transitions as shown in 
our theoretical results. 

To further support our analysis, we perform double¬ 
derivative analysis of critical points of STO dielectric 
function at 300 K as representatives. The structures in 
the optical spectra (e) can be resolved into multiple crit¬ 
ical points using standard analytical line-shapes [26, 36- 
38] given by 

e{u)=C-Ae i *(u-E + tT) n (9) 

Here the critical point parameters A, E, T and rep¬ 
resent amplitude, threshold energy, broadening and exci¬ 
tonic phase angle respectively. The exponent n depends 
on the dimension of the critical point in consideration. 

Double derivative fit of the (pseudo-) dielectric func¬ 
tion of STO and SrTq _ r Nb x () 3 , 0<a;<0.005 have been 
performed for all the measured temperatures. The equa¬ 
tions used for fitting are derived from Equation 9 and are 
given below. 

For a discrete excitonic critical point lineshape 

= n(n - 1 )Ae**(u -E + iT) n ~ 2 , n = -1 (10) 

and for a two dimensional (2D) critical point lineshape 

44 =Ae i *{u J ~E + iT)- 2 (11) 

au> z 

A representative example of such fitting analysis for 
the critical points of STO at 300 K is presented here 
in Fig. 5. The fitting parameters for the best fit (with 
MSE of 4.287) are given in Table I. The numbers in the 
parentheses indicate 95% confidence limit. 

The evolution of the optical spectra of SrTii-^Nb^Oa 
with different x (0.0001, 0.0005, 0.001 and 0.005) is 
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FIG. 7. (Color online) Experimental critical points and 
comparison with theory for SrTii-^Nb^Oa, x = 0.005. (a) 
(ei) and (£ 2 ) for the lowest (4.2 K) and highest (350 K) mea¬ 
sured temperatures, (b), (c) Temperature dependent differ¬ 
ences 5t(£z) = (£i)(T) — ( £i)(T = 350 K),i = 1,2. (d) Plots of 
£1 and £2 calculated using GW-BSE for SrTii-^NbajOs, x = 
0.005 (all graphs are blue-shifted by 53 meV to match exper¬ 
imental data). The green vertical dashed line are determined 
using second derivative analysis. 


shown in Fig. 6. The room temperature plots of (si) and 
(£ 2 ) (Fig. 6(a)) exhibit a strong Drude tail for the sam¬ 
ple with highest doping (x = 0.005). Differences 5d{ei) 
and Sd{s 2 ) between the PDF of STO and SrU^^Nb^Oa 
(S d (£i) = {£i){x) - ( £i)(x = 0 ),i = 1,2) highlight the 
doping-dependent evolutions as shown in Figs. 6(b) and 
6(c), respectively. The changes are manifested as a Drude 
tail and several negative peaks in 5d(e 2 )- The peak po¬ 
sitions gradually shift to higher energies with increasing 
doping. 

Another important observation is the evolution of the 
excitonic peaks and the corresponding rearrangement of 
the spectral weight with the variation of x (Fig. 6(c)). 
With the increase of x from 0 to 0.001 it is seen that the 
spectral weight is redistributed only between the energy 
range from ~3 to 6.5 eV. For example, the negative peak 
at ^3.9 eV is compensated by the positive contributions 
from above ~4.8 eV. However, in the case of x = 0.005 
the onset of the Drude tail is seen along with a change 
in trend of the spectral weight transfer. Spectral weight 
redistribution is not confined to the energy range from ~3 
to 6.5 eV anymore as before, conspicuously changing to 
negative above ~ 4.8 eV. This could be an indication that 
SrTii-^Nb^Oa is still in the insulating state for x < 0.001 
despite the presence of doping. 

In Figure 7(a) the PDFs at 4.2 K and 350 K are plot¬ 
ted for x = 0.005. The PDF temperature difference 
plots in Figs. 7(b) and 7(c) (St{s 1 ) and St{s 2 ), respec¬ 
tively) show signatures of a Drude tail and sharper indi¬ 
rect band-edge features. The enhancement of the indi¬ 
rect band-edge feature could be due to stronger electron- 
phonon coupling [39]. Similarly to the case of STO, ex¬ 
citonic peaks are seen as prominent features in St{s 2 )- 
From double derivative analysis these peaks are found 
to be at 3.75, 4.76 and 6.11 eV as compared to 3.75, 
4.67 and 6.11 eV, respectively, in case of STO. The blue- 
shift in the case of the localized exciton Ex 2 , together 
with broadening and decrease of heights of all the peaks, 
can be interpreted as a reduction of the excitonic effect 
due to the presence of Nb-dopant free charge carriers 
which increases the screening and reduces the e-h in¬ 
teraction. Systematic measurements performed also on 
SrTii_ 3 ,Nb;j 03 , x = 0.0005 (which is an intermediate 
value) show similar trend. This mechanism is confirmed 
by GW-BSE simulations (Fig. 7(d)) with doping and 
these results represent further evidence of the dominant 
excitonic nature of these peaks. The dielectric function 
calculated from GW-BSE for SrTii-^Nb^Oa, x = 0.005 
are plotted in Fig. 7(d). There is very good agreement 
with the experimental PDF results in terms of excitonic 
peak positions as well as other prominent features. 

IV. CONCLUSION 

In summary, the evidence of anomalous excitonic ef¬ 
fects in the optical spectra of STO in the form of a 










Wannier -like bound exciton at 3.75 eV, resonant excitons 
at 4.67 eV and 6.11 eV at 4.2 K, as well as their evolu¬ 
tion with temperature, doping-modulated screening and 
correlation, are reported. Remarkably the Wannier- like 
bound exciton is more delocalized than the resonant exci¬ 
ton. Theoretical calculations further show the indispens¬ 
able role e-h and e-e interactions, yielding to excitonic 
effects, play in understanding the dielectric function as a 
whole. These observations of prominent and novel corre¬ 
lated phenomena in STO may pave the way for correct 
understanding of not only STO per se but also mate¬ 
rial systems where it is an integral component such as 
heterostructures. 
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